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We report on a search for bottom squarks (b) produced in pp collisions at ^/s — 1.8 TeV using the 
D0 detector at Fermilab. Bottom squarks are assumed to be produced in pairs and to decay to the 
lightest supersymmetric particle (LSP) and a b quark with a branching fraction of 100%. The LSP 
is assumed to be^the lightest neutralino and stable. We set limits on the production cross section 
as a function of b mass and LSP mass. 
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Supersymmetry (SUSY) is a hypothetical fundamen- 
tal space-time symmetry relating bosons and fermions 
0. Supersymmetric extensions to the standard model 
(SM) feature as yet undiscovered supersymmetric part- 
ners for every SM particle. The scalar quarks (squarks) 
5l and qR are the partners of the left-handed and right- 
handed quarks, respectively. These are weak eigen- 
states, and can mix to form the mass eigenstates, with 
qi = §z,cos0 + qnsmO for the lighter squark, and the or- 
thogonal combination for the heavier squark §2 ■ In most 
SUSY models, the masses of the squarks are approxi- 
mately degenerate. But in some models, the lighter top 
and bottom squarks could have a lower mass than the 
other squarks because of the high mass values of the top 
and bottom quarks. In particular, lighter bottom squarks 
could arise for large values of tan/3, the ratio of the vac- 
uum expectation values of the two Higgs fields in the 
minimal supersymmetric standard model. 

We report the results of a mixing-independent search 
for bottom squarks produced in pp collisions at *Js =1.8 
TeV. Squarks are produced in pairs by QCD processes 
with the production cross section depending on the mass 
of the squark but not on the mixing angle 9. We search 
for events where both squarks decay to the lightest neu- 
tralino x? via b — ► X? + ^ and assume that the Xi 1S 
the lightest supersymmetric particle (LSP) and stable. 
This should be the dominant decay channel provided that 
the mass of the squark (mg) is larger than the combined 
masses of the b quark and LSP (mLSp); therefore we as- 
sume its branching fraction is 100%. This yields a final 
state consisting of two b quarks and two unobserved sta- 
ble particles resulting in missing transverse energy ($ T ) 
in the detector. In this paper, we give limits on the 
squark pair production cross section for different values 
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. 1. The expected distributions of If T for mj values of 
and 100 (b) GeV/c 2 , for the indicated values of ttilsp 



jets topology. This channel was previously used to set 
limits on the mass of the top squark, which was assumed 
to decay t — > Xi + c Pi • The other three channels in addi- 
tion required that at least one jet has an associated muon, 
thereby tagging b quark decay, and were used to set limits 
on a charge 1/3 third generation leptoquark for the de- 
cay LQ — > v T + b H . We use identical data samples and 
event selections for the bottom squark limits presented 
in this paper. For all channels, the presence of signif- 



of mi and m LS p. Limits on the cross section are used to i can t $ T is used to identify the non- interacting LSPs. 
exclude a region in the (mLSP, mj) plane. Limits [g| from 



the CERN e + e~ collider (LEP) experiments depend on 
the Z/7-to-squark coupling, which is a function of the 
mixing angle. For maximal coupling, the LEP exclusion 
region can extend to the kinematic maximum; for exam- 
ple, to about 85 GeV/c 2 at ^ = 183 GeV. 

The data used for our analysis were collected during 
1992-1996 by the D0 detector g] at the Fermilab Teva- 
tron Collider. The D0 detector is composed of three ma- 
jor systems: an inner detector for tracking charged par- 
ticles, a uranium/liquid argon calorimeter for measuring 
electromagnetic and hadronic energies, and a muon spec- 
trometer consisting of a magnetized iron toroid and three 
layers of drift tubes. The detector measures jets with 
an energy resolution of approximately a/E = 0.8/yE 
[E in GeV) and muons with a momentum resolution of 
2 



distribution for two val- 
b and different ttilsp 0- Our requirement that 
— 40 GeV reduces the acceptance for small val- 



<r/P = [( al8( p P ~ 2) ) 2 + (0.003P) 2 ] 1 / 2 in GeV/c). $ T is 
determined by summing the calorimeter and muon trans- 
verse energies, and is measured with a resolution of a = 
1.08 GeV + 0.019-(£|£r|) §. 

Four channels are combined to set limits on the pro- 
duction of bottom squarks. The first required a $ T and 



Figure 1 shows the expected $ 
ues of m 
E T > 35 

ues of the mass difference m^-rriLSP ■ Backgrounds arise 
from events where neutrinos produce significant for 
example, in lU+jets events, where W — > Iv. 

Events for the ^ T +jets channel were collected using a 
trigger that required $t> 35 GeV. The offline analysis 
required two jets (E^ > 30 GeV), $ T > 40 GeV, and 
no isolated electrons or muons. Events had to have only 
one primary vertex to assure an unambiguous calcula- 
tion of $ T . To eliminate QCD backgrounds, additional 
cuts were made on the angles between the two jets, and 
between jets and the direction of the ]fi T . Data with an 
integrated luminosity of 7.4 pb" 1 , satisfying the above 
selection criteria, yielded three candidate events. Back- 
ground was estimated to be 3.5±1.2 events, with 3.0±0.9 
events from W boson decays and 0.5 ± 0.3 events from Z 
boson decays S. 

The trigger for the muon channels required either two 
1ow-pt muons > 3.0 GeV/c), or a single low-pr muon 
and a jet with Et > 10 GeV, or a high-py muon (pt^ > 15 
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GeV/c) and a jet with Et > 15 GeV. Integrated lumi- 
nosities of 60.1 pb -1 , 19.5 pb -1 , and 92.4 pb -1 respec- 
tively were collected using the three muon triggers. The 
offline analysis used muons in the pseudorapidity range 
\rjnl < 1.0 and p T > 3.5 GeV/c, while jets were required 
to have Et > 10 GeV. For events with two muons, each 
muon had to be associated with its own jet. In single 
muon events, the muon was required to be associated 
with a jet, and an additional jet with Et > 25 GeV was 
also required. To remove QCD backgrounds, events were 
selected with Jp T > 35 GeV and an azimuthal angular 
separation between the $ T and the nearest jet of > 0.7 
radians. For the single muon channels, backgrounds from 
W boson decays were reduced by cuts on muon-jet corre- 
lations, while background from top quark production was 
minimized by cuts on the scalar sum of jet Et- After im- 
position of all selection criteria, two events remained in 
the data. 

We considered background contributions to the muon 
channels from tt and W and Z boson decays Q. Top 
quark events have multiple b quarks and $ T , and we es- 
timated that 1.4 ± 0.5 tt events remained in our sample. 
W and Z events have ~$ T from W — > Iv or Z — > vv. 
They can also have muons near jets that can mimic b 
quark decays when a prompt muon overlaps a jet, or a 
jet fragments into a muon via a c quark or a tt/K de- 
cay. We estimated there were 1.0 ± 0.4 W boson events 
and 0.1 ± 0.1 Z boson events in the sample. The total 
background for the muon channels was therefore 2.5±0.6 
events. 

Combining the four channels yields five events, with 
a total estimated background of 6.0 ± 1.3 events. We 
set limits on the cross section by combining the detec- 
tion efficiencies and integrated luminosities for the dif- 
ferent channels. We calculate the detection efficiency us- 
ing Monte Carlo (MC) generated acceptances 0, multi- 
plied by trigger and reconstruction efficiencies obtained 
from data [5,6]. The total efficiencies for different squark 
and neutralino masses are summarized in Table I. Us- 
ing a muon to tag b quark decays reduces the efficiency 
for those channels, but their higher integrated luminosi- 
ties yield a sensitivity comparable to that of the $ T +jets 
channel. Including systematic errors and statistics for 
the MC, the total uncertainty on the combined efficiency 
varies between 8.6% and 29%, depending on the assumed 
masses. The jet energy scale dominates the systematic 
error for m~ b = 70 GeV/c 2 , while uncertainties on the 
muon trigger and reconstruction efficiency dominate at 
higher squark masses. The 95% confidence level (C.L.) 
upper limits on the pair production cross section are de- 
termined using Bayesian methods, and include the sys- 
tematic uncertainty on the efficiency and a 5.3% uncer- 
tainty in the integrated luminosity. The resulting upper 
limits are given in Table I for different values of m b and 

We use the program PROSPINO [|| to calculate the bot- 
tom squark pair production cross section as a function of 
m b . The cross section is evaluated assuming a renor- 



TABLE I. Total efficiencies for different and tulsp val- 
ues for the four channels, and 95% C.L. limits on the produc- 
tion cross section obtained by combining all channels. 
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FIG. 2. The 95% C.L. exclusion contour in the (mLSP,^) 
plane. Also shown are the results from the ALEPH experi- 
ment at LEP for minimal (9 — 68°) and maximal (9 — 0°) 
coupling [2]. 

malization scale fi = m~. The program includes next- 
to-leading order diagrams, and uses CTEQ4M parton dis- 
tribution functions fa]. For any given mg, we determine 
the value of mLSP where our 95% C.L. limit intersects 
the theoretical cross section. The excluded region in the 
(wLSPj TO h) pl ane is shown in Fig. 2. We exclude values 
of m~ b below 115 GeV/c 2 for to lsp < 20 GeV/c 2 . For 
rn b = 85 GeV/c 2 , we exclude the region with mLSP < 47 
GeV/c 2 . Also shown are limits (2| from ALEPH for 
y/s = 181 — 184 GeV. For most allowable values of tulsp, 
they exclude the region with m b < 83 GeV/c 2 , assuming 
maximal coupling (9 = 0°) |i~0| ]. 

In conclusion, we observe five candidate events con- 
sistent with the final state bb+$ T - We estimate that 
6.0 ± 1.3 events are expected from tt and W and Z bo- 
son production, and find no excess of events that can be 
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attributed to bottom squark production. We interpret 
our result as an excluded region in the (ttilsPj'tI;;) plane. 
This result is independent of the mixing between hi, and 
~b R . 
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